Electronic beam shifts in monolayer graphene superlattice 
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Electronic analogue of generalized Goos-Hanchen shifts is investigated in the monolayer graphene 
superlattice with one-dimensional periodic potentials of square barriers. It is found that the lateral 
shifts for the electron beam transmitted through the monolayer graphene superlattice can be negative 
as well as positive near the band edges of zero-fc gap, which are different from those near the band 
edges of Bragg gap. These negative and positive beam shifts have close relation to the Dirac point. 
When the condition qAdA = —qadB = rmr (m = 1,2,3...) is satisfied, the beam shifts can be 
controlled from negative to positive when the incident energy is above the Dirac point, and vice 
versa. In addition, the beam shifts can be greatly enhanced by the defect mode inside the zero-fc 
gap. These intriguing phenomena can be verified in a relatively simple optical setup, and have 
potential applications in the graphene-based electron wave devices. 
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I. INTRODUCTION 

Since the first synthesis of graphene in 2004 [H-Q , elec- 
tronic analogies of optical phenomena ffll^ , such as re- 
flection, refraction, polarization, interference, waveguid- 
ing, focusing and cohimation have inspired many inter- 
esting and new concepts on "Dirac electron wave optics" 
to design the graphene-based nanoelectronic devices, due 
to the fink between the Klein tunneling |15| and negative 
refraction [l^. The close relation between electron and 
light results from the quantum-mechanical wave nature 
of electrons. Among all these optical-like phenomena, 
Goos-Hanchen shift, referring to the lateral shift for to- 
tally reflected light beam at a single dielectric interface 
[l3| , has been extended to the field of electronics [l^ and 
currently investigated in various graphene nanostructures 
[19I - I23I . The electronic lateral shifts result in a remark- 
able ie^ /h conductance plateau in p-n-p graphene junc- 
tion [l9| , and further provide an alternative way to realize 
valley beam splitter in strained graphene (22l . [23j . 

Motivated by the experimental progress on the 
construction of graphene superlattice, the electronic 
bandgap and transport in graphe ne superlattice have 
attracted much attention |2J-|31[. Interestingly, the 
graphene superlattice is analogous to one-dimensional 
photonic crystal containing positive-index and negative- 
index materials jsO, 31]. thus there exists zero-averaged 
wave number (zero-A:) bandgap in graphene superlat- 
tice. Like the zero-n gap in the photonic crystals with 
negative-index materials [s^ . this so-called zero-fc gap, 
which is associated with the Dirac point in graphene su- 
perlattice m, im , is difi'erent from Bragg gap and thus 
robust against the lattice parameters and incidence an- 
gles. Besides the unique electronic band gap and trans- 
port in the graphene superlattice, the beam supercolli- 
matation is also achieved in two-dimensional graphene 
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superlattice, utilizing the highly developed concepts in 
optics 0. 

In this work, we will investigate the generalized Goos- 
Hanchen shifts of electron beam in transmission through 
monolayer graphene superlattice with one-dimensional 
periodic potentials of square barriers. The beam shifts 
discussed here has nothing to do with the evanescent 
wave and total rcfiection, so we term it as generalized 
Goos-Hanchen shift, since the lateral shift, resulting from 
the finite width of the beam, keeps the main features of 
the Goos-Hanchen shift. The behaviors of beam shifts 
are quite different near the band edges of the zero- A: gap 
and Bragg gap. What we emphasize here is that the un- 
usual feature of beam shifts around the Dirac point in 
the graphene superlattice, that is to say, when the in- 
cident energy is below the Dirac point, the beam shifts 
can be positive and vice versa. Moreover, there will ex- 
ists transition between the negative and positive beam 
shifts, depending on the condition qacIa ~ —qsdB = "itt 
(to = 1,2,3...). When the defect in the graphene super- 
lattice is further induced, the beam shifts inside the zero- 
fc bandgap can be greatly enhanced near the defect mode, 
whose location depends on the incidence angles at a fixed 
incident energy. Furthermore, the beam shifts can also 
be controlled by adjusting the potential height of defect 
in the graphene superlattice. Consequently, the tunable 
negative and positive beam shifts will provide the poten- 
tial applications in graphene-based electron wave devices 
(e.g., beam switch, beam splitter and spatial modulator). 



II. THEORETICAL MODEL 



Consider a graphene-based superlattice with one- 
dimensional periodic squared potential barriers. The 
schematic structure for the monolayer graphene super- 
lattice is shown in Fig. [TJ where the parameters are the 
potential barrier height Va and width (Ia, the potential 
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FIG. 1. (Color online) (a) Schematic diagram for beam shift 
of the transmitted electrons through a graphene superlattice 
with defect C. (b) The profiles of the potentials. 



well height Vb and width ds, and the defect is denoted 
by the potential height Vc and width dc- In what fol- 
lows, we will first study the lateral shift in the graphene 
superlattice (AB)^ , and discuss further the influence of 
defect mode in the structure of {AB)^C{BA)^ . 

In general, the Hamiltonian of electrons inside a mono- 
layer graphene, in the vicinity of the K point and in the 
presence of a potential V{x), is given by the following 
massless Dirac-like equation. 



H = -ihvFa-V + V{x), 



(1) 



where the Fermi velocity vp ~ lO^m/s, and t? = (ax, cry) 
are the Pauli matrices. Due to the translation invariance 
in the y direction, the solution of above equation for a 
given incident energy E and potential barrier Vj can be 
presented as ^(a;,y) = 'i>{x)e^''^y with 
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iqj3 
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where the spinor wave function ^ = 



(2) 
has the 



pseudospin component 4*+ and on the A and B sub- 
lattice, kj = (_B — Vj)/hvF, ky and qj arc the y and x 
components of wave number, qj = sign(fcj)(fc| — kyY^^ 
for kj > ky, otherwise qj ~ i{ky — kj)^^'^, and aj (bj) 
is the amplitude of the forward (backward) propagating 
wave. Noting that when E < Vj, the wave vector is op- 
posite to the direction of electron's velocity. This results 
in the negative refraction in p-n junction, and Veselago 
lens which is relevant to the negative Goos-Hanchen 
shifts in graphene barrier plj . The wave functions at any 
two positions x and x + Ax inside the jth potential can 



be related by the transfer matrix [3ll |: 



cos((jj Ax — 6j ) ■ sin{qj Ax) 
cos Oj cos Oj 

■ s\n{qj A.x) cos{qj Ax-^Oj ) 



(3) 



with 6j ~ SLTCsm{ky/kj). After lengthy but direct cal- 
culations, the reflection and transmission coefficients are 
found to be 



r{k 

t{ky 



(X22e-*^o + xiie*^*) - xi2e^(^*-''«) - X21 
2 cos 6*0 



(X22e- 



Xiie*^*) - xi2e' 



X21 



,(4) 
,(5) 



where and Ot are incidence and exit angles (see Fig. [T|) , 
and Xij(i,j = 1,2) is the matrix element of total transfer 

matrix, Xn ~ YijLi connecting the incident and exit 
ends, and N is the total number of the layers of graphene 
superlattice. 

To calculate the generalized Goos-Hanchen shift, we 
consider an electron beam, 
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dkyA{ky)e'^''^''+''yy^ 



1 



(6) 



where k^ and ky are the x and y components of wave 
number in the incident region, k^ = (fc^ — kyY^^, 9 = 
arcsin (fey/A:) and k = E/hvp with assuming V = 0. The 
angular-spectrum distribution is assumed to be, but not 
necessarily, A{ky) = Wy exp [^{w'^/2){ky — kyof'] around 
kyQ corresponding to the incidence angle 6q, where Wy = 
w sec 60, and w is the half beam width at waist. When 
the electron beam is well-collimated, there will be a nar- 
row distribution of ky values around kyo. Thus, we can 
expand k^ and 9 to first order around kyo and sub- 
stitute into Eq. ([6]). By evaluating the Gaussian in- 
tegral, we finally obtain the spatial profile of the in- 
cident beam at x ~ 0, which leads to the two com- 
ponents as follows, = exp [— (y — j/™)^/2Wy], with 
yf = and = ~0'{kyo). Noting that 0'{kyo) means 
the values of derivative with respect to ky taken at 
ky = kyo- Obviously, the separation of the two centers is 
<^o = |y™ — y™| = —S'ikyo). For the convenience, we can 
consider the potential in the transmitted region is the 
same as that in the incident one, so that the transmitted 
beam can be then expressed as 
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V2^ J- 



dkyt(^ky) A(^ky )c 



i[k:c{x-dT}+kyy] 



1 



(7) 

where dp is the total length of the graphene superlat- 
tice and the complex transmission coefficient t(ky) for 
the relevant plane wave component can be expressed 
in term of amplitude p{ky) and phase ip{ky), t{ky) = 
p{ky) exp [iip{ky)]. As a matter of fact, when the station- 
ary phase approximation is applied, it is assumed that 
the phase ip{ky) of the transmission coefficient t{ky) is 
linearly dependent on ky and the amplitude of the trans- 
mission keeps almost constant, thus it is reasonable to 
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III. RESULTS AND DISCUSSIONS 
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FIG. 2. (Color online) (a) Dependence of beam shifts and 
transmission probability for transmitted electron beams in 
graphene superlattice (AB)^^ on the incident energy E, where 
00 = 20°, Va = 50 meV, Vb = meV, ds/dA = 1, = 20 
nm (solid blue line), and dA = 30 nm (dashed red line); (b) 
and (c) are the corresponding electronic band structures for 
dA = 20 nm and dA = 30 nm. The dot lines denote the 
location of the Dirac point Ed = 25 meV, and the dashed 
lines present the "light cones" of the incident electrons. 



expand the f{ky) in Taylor series at fcj^Oj smd retain up to 
the first-order term. Similarly, kx and 9 can be expanded 
to first order around fcyo- After substituting all into Eq. 
([7]), we can calculate the Gaussian integral, and obtain 
the two components of the transmitted beam at x = dx, 
exp[-(y-y±)V2w^], with = -v'{kyo) and 
~Lp'{kyQ) — 9'{kyo). Noting that the separation of 
the two centers is also Sq, which is the same as that for 
incident beam. 



Comparison of the incident and transmitted beams, 
the displacement of two components are the same, a± = 
_ _ gQ the lateral shift of the 

electron beam transmitted through the graphene super- 
lattice can be defined as [l^, A = ((t+ +(t^)/2, and thus 
is given by [H [U 



A = V(^yo) = -Im(^lnt), (8) 



where d/dkyo denotes the derivative with respect to ky 
evaluated at ky = kyo. 



Fig. [2] (a) shows the typical beam shifts in transmis- 
sion through the graphene superlattice (AS)^°, where 
6*0 = 20°, Va = 50 meV and Vb = meV, ds/dA = 1, 
dA = 20 nm (solid blue line), and Fig. [2](b) and (c) illus- 
trates the corresponding electronic band-crossing struc- 
tures in which Ed = 25 meV. Surprisingly, the beam 
shifts near the band edges of the zero-fc gap and Bragg 
gap are quite different. As shown in Fig. [2l the beam 
shifts can be negative as well as positive near the band 
edges of the zero-A: gap, which depends on the incident 
energy below or above the Dirac point, Whereas, 
the beam shifts near the band edges of the Bragg gap are 
positive. Remarkably, the enhanced beam shifts in this 
case are sensitive to the lattice constants, dA and ds, 
because the Bragg gap depends strongly on the lattice 
constants. 

To understand beam shifts in such a graphene super- 
lattice in term of electric band gap, we write down the 
electronic dispersion at any incidence angles for the infi- 
nite superlattice N oo, based on the Bloch's theorem, 
in the following form: 

/n ^ ^ I i , \ COS (6 A ^ Ob) - 1 

COS [Pxh) = cos [qAdA + qBdB) H 7 7 

cos Oa cos Ob 

X sin [qAdA) sin (qBdB), (9) 

where /3x is the x component of Bloch wave vector, and 
A = dA + dB is the length of the unit cell. When qAdA = 
—iBds "iTT-^ ("^ = lj2, 3...), and Oa 7^ 0, there is no 
real solution for thus there exists zero-fc bandgap 31 



The Dirac point, the location of the touching point of 
the band, is given by q^d^ = — te^s at Oa = 0. In 
the simple case of 7^ and Vg = 0, the Dirac point is 
located at = Va/(1+c'b/'^a)- For instance, the Dirac 
point is calculated as Ed = 25 meV with the parameters 
Va = 50 meV and dB/dA ~ 1. From the electronic band 
structures in Fig. [2](b) and (c), one can predict that the 
group velocity is always negative (positive), when the 
incident energy is below (above) the Dirac point, based 
on the electronic dispersion relation. So this provides the 
intuitive explanation of the negative and positive beam 
shifts in this case. 

More interestingly, according to the electronic disper- 
sion (O, when the condition 



qAdA = -qsdB = mn, (?7i = 1, 2, 3. 



(10) 



is satisfied, the zero-fc gap will be closed, and a pair of 
new zero-fc states emerges from ky = 0. Actually, the 
extra Dirac points, located at the incline angles, for in- 
stance ±fcy 7^ 0, are only dependent on the lattice con- 
stant A = dA + dB and the ratio of dB/dA [3l|. This 
condition PH)) is the same as that in one-dimensional 
photonic crystals containing negative-index and positive- 
index materials (33| . Comparison of the electronic band 
structures for different lattice constants demonstrates 
that the zero-fc gap can be close and open periodically 
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FIG. 3. (Color online) Dependence of the beam shifts on the 
incidence angle near the extra Dirac point, where (a) E > 
Eu and (b) E < Eq- (c) and (d) are the corresponding 
phase shifts, (e) Electronic band structure for the graphene 
superlattice with (Ib ~ dA = 120 nm. 



with the increasing the lattice constants. In the follow- 
ing discussion, we shall focus on the beam shifts near the 
band-crossing structure, when the graphene superlattice 
satisfies the condition ([TO)). which offers the freedom to 
control the lateral shifts around the Dirac points at non- 
normal incidence angles. 

Fig. [3] (a) illustrates the dependence of the beam shifts 
on the incidence angle in the graphene superlattice, when 
the condition PU)) with m = 1. Above the Dirac point 
{E > Ejj), there exists transition from the negative beam 
shifts to the positive ones with increasing the incidence 
angle ^o- On the contrary, the beam shifts changes from 
positive to negative below the Dirac point {E < Eu). To 
confirm these intriguing phenomena, the corresponding 
phase shifts of the transmitted beam are depicted in Fig. 
[3](c) and (d), in which the slopes of the phase shifts sug- 
gest the negative and positive lateral shifts. As a matter 
of fact; the beam shift can be understood from the re- 
shape of the transmitted beam whose plane components 
undergo the different phase shifts. These beam shifts 
are quite different from the previous ones, for example, 
d-A = d-B = 20 nm. In that case, QAdA = —qBds < ti", so 
the beam shifts near the band edges of the zero-fc gap are 
always negative and positive, depending on whether the 
incident energy is above or below the Dirac point. How- 
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FIG. 4. (Color online) Effect of the defect on beam 
shift and transmission property in graphene superlattice 
(AB)"'C(B4)"', where dA = dg = 20 nm, Va = 50 meV, 
Vs = meV, dc = 70 nm, (a) Vc = 50 meV, fo = 15° 
(solid blue line), 6*0 = 20° (dashed red line) Go = 25° (dotted 
black line), (b) Bq = 20°, Vc = 40 meV, (solid black line), 
Vc = —40 meV (dotted red line), Vc = 45 meV (dashed blue 
line), and Vc ~ —45 meV (dash-dotted green line). 



ever, there is extra Dirac point at ky ^ 0, as shown in 
Fig. [3](e). Due to the periodical appearance of the Dirac 
point, the group velocity can be negative and positive 
with increasing the incidence angle. As a consequence, 
the beam shifts can be controlled from negative to posi- 
tive, and vice versa. 

Next, we shall investigate the influence of the defect 
mode on the beam shift. It is found in Fig. |l](a) that the 
beam shift at the incident energy near the defect mode 
will be greatly enhanced at different incidence angles, and 
there exists defect mode with its incident energy inside 
the zcro-fc bandgap. Actually, the beam shift is almost 
zero in this band gap, except the near the corresponding 
energy of defect mode. Fig. |4] (b) further demonstrates 
that the enhanced lateral shift can be controlled from 
negative and positive by adjusting the potential of de- 
fect, Vc- It should be noted that the magnitude of the 
enhancement near the defect mode can be much larger 
than that near the edges of band gap. Therefore, the 
control of beam shifts over the defect provides more flex- 
ibility. 
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IV. CONCLUSION 



To summarize, we have investigated the generahzcd 
Goos-Hanchen shifts for the electron beam transmitted 
through the monolayer graphene superlattice. The beam 
shifts near the edges of zero-fc gap can be negative as well 
as positive, which depends on the incident energy above 
or below the Dirac point. When qAdA = —QBds = ^titt 
(m = 1,2,3...) is satisfied, there can exist transition be- 
tween the negative and positive beam shifts. In addi- 
tion, the negative or positive beam shifts can be enhanced 
and controlled by the defect mode when the defect is in- 
duced in the graphene superlattice. These beam shifts in 
the graphene superlattice are applicable to design various 
graphene-based electron wave devices. 

Last but not least, recent progress on analogous phe- 
nomena related to Dirac point in photonic crystals 
[33I [s^ l or negative-zero-positive index metamaterial (35| 
suggest that many exotic phenomena in graphene can 



be simulated with relatively simple optical benchtop ex- 
periments. Therefore, we also hope that the electronic 
beam shifts in graphene superlattice presented here can 
be tested in the analogy of the photonic crystals contain- 
ing negative-index and positive-index materials [s^ . 
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